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ABSTRACT 


In this paper, two optimization methods are used to adjust the gain values for 
the cascade PID controller. These algorithms are the butterfly optimization 
algorithm (BOA), which is a modern method based on tracking the 
movement of butterflies to the scent of a fragrance to reach the best position 
and the second method is particle swarm optimization (PSO). The PID 
controllers in this system are used to control the position, velocity, and 
current of a permanent magnet DC motor (PMDC) with an accurate tracking 
trajectory to reach the desired position. The simulation results using the 
MATLAB environment showed that the butterfly optimization algorithm is 
better than the particle swarming optimization (PSO) in terms of 
performance and overshoot or any deviation in tracking the path to reach the 
desired position. While an overshoot of 2.557% was observed when using the 
PSO algorithm, and a position deviation of 7.82 degrees was observed from 
the reference position. 


This is an open access article under the CC BY-SA license. 


Corresponding Author: 


Kareem Ghazi Abdulhussein 

Department of Electrical Power Engineering Techniques 
Middle Technical University 

Babylon, Iraq 

Email: engkarim1984@ gmail.com 


1. INTRODUCTION 

A permanent magnet DC motor (PMDC) is a simple structure type of DC motor in which the field 
windings have been replaced by permanent magnets and have many applications so It operates as the same 
basic principle as a shunt connected motor but the difference is that the permanent magnet generates the 
required flux instead of field windings [1], [2]. 

Permanent magnet DC motors provide less output power than DC shunt motor because the flux 
generated from permanent magnets is less than what can be achieved with field windings so most types of 
small DC motors are permanent magnet DC motors and these motors usually operate at high velocity and 
little torque [1], [3]. This motor has many applications such as its use in car front wipers and moving 
windows for cars and household appliances such as food mixers and others [4], [5] .1n addition to the 
importance of using it in CNC machines, robotics, and electric vehicles [6], [7]. 

The main problem statement in this paper is how to get the best parameter values for cascade PID 
controller which gives the best accurate results for tracking trajectory the reference position and reaching the 
desired position at a regular velocity. To address this problem statement, a comparison was made between the 
BOA and PSO methods for extracting the parameters of the cascade PID controller, in addition to using one 
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of the performance criteria which is integral time absolute error (ITAE) to reduce the error between input and 
output for this system. 

The cascade PID controller which is used in this paper consists of three controllers, the current 
controller is designed as an inner loop, and both position and velocity are outer controllers [8]. P, PI, and PI 
controllers are used instead of a PID controller for position, velocity, and current respectively [9]. 

The purpose of using the cascade control unit is due to several reasons, the most important of which 
is to reduce or reject disturbance and return to the steady-state [10]-[12]. There are several methods to tuning 
the parameters of PID and use it in PMDC motor such as Ziegler-Nichols Method (Z-N), Cohen-Coon 
method, artificial neural, network fuzzy logic [13], [14], particle swarm optimization (PSO) [15], [16], 
genetic algorithm [17], and butterfly optimization algorithm (BOA) [18]. 

To compare this work with similar works, a literature review was made. In 2014 Md Mustafa et al 
introduced a BLDC speed control system using GA [19]. In 2015, Taha et. al. used three methods to control 
the cascade control system [9]. In 2018, Wisam et. al. introduced a system for controlling PMDC velocity 
using the GA and DS methods [20]. in 2019, Fadhil et. al. used a fractional PID controller to control PMDC 
speed based on PSO [21]. In 2021, Ahmed et al introduced a system to control the speed and position of the 
servo motor [22]. 

This paper is organized as follows: In the second section, the mathematical model of the PMDC 
motor is represented and explain the general structure of the system, in the third section, the two tuning 
methods are explained which are PSO and BOA in addition to the use of ITAE. The fourth section includes 
the results and comparison, while the final section includes a conclusion. 


2. MATHEMATICAL MODEL OF PMDC MOTOR AND GENERAL STRUCTURE 
2.1. Mathematical Model of PMDC Motor 

The equivalent circuit for a PMDC motor shown in Figure 1 consists of an armature circuit 
consisting of inductance (La) and resistance (Ra) connected on series. An electromotive force (Ea) is 
generated as a result of cutting the lines of flux generated by the permanent magnet, and it has a direction 
opposite to the direction of the source voltage. The mechanical part consists of moment inertia (Jm) and 
friction coefficient (Bm). The other parameters are torque constant (K+) and back emf constant (K,). The block 
diagram of PMDC motor with the position output showed in Figure 2. 


Vals) 





Figure 2. Block diagram of permanent-magnet DC motor [23] 
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From the block diagram, The Mathematical Model of PMDC motor with position output can be 
listed in the (1) to (4) [23]. 


Va (©) = ea (t) + Raia (©) + La Sia (t) (1) 
ea (t) = Ke W(t) (2) 
Tm Ct) — Th = Jim <= @m(t) + Bam (t) (3) 
Ta (© = Kiia (© (4) 


By using Laplace transformation for the (1) to (4) we obtain. 


Va (s) = Ea (s) + Ra la(s) + SLala(s) (5) 
Ea(s) = Ke W@m(S) (6) 
Tm (S) — TL = SJm@m(s) + Bn @m(S) (7) 
Tm (S) = Ky la(s) (8) 


The overall transfer functions are defined for velocity and position control of the PMDC motor, respectively. 


om(S) _ ke = (9) 
Va (s) J LaS2+(JRa +BLa ) S+BRa + ke” 





o(s) — ke 
Va (s)  JLaS3 +(JRa +BLa ) S2 +( BRa + ke” )S 





(10) 


Where ke equal to K, [24]. 
Table 1 lists the main components and their values of PMDC motor Parameters. 


Table 1. The PMDC Parameters [9] 


Motor Parameters Value 

Torquen constant K,= 2.35N.m/A 
Armature inductance L, = 2.61*10°H 
Armature resistance R,2.61Q 

Inertia constant Jm = 0.068kg.m? 
Friction coefficient Ba = 0.008N.m.s/rad 
Back emf constant K, = 2.35V.s/rad 
Load torque T= 17.6N.m 
Applied voltage Va = 230v 


2.2. General structure of system 

Figure 3 represents the general structure of the system, which consists of three controllers for 
current, velocity, and position. The output from the position controller represents the reference velocity, the 
output from the velocity controller represents the reference current, while the output from the current 
controller represents the control voltage (vc) and it is compared with a triangular waveform signal that is 
generated internally in the PWM circuit, producing a KPWM gain that simply represents the PWM controller 
and the dc-dc switching mode converter (1.e. H bridge). KPWM can be represented by the following equation 
in Laplace transform after making linearized of the dc-dc converter [24]. 


Va(s) = KPWM = V,(s) (11) 
One of the benefits of the cascade control system, in addition to the benefits mentioned previously, is 
the ability to put limits on the reference signals for the system in order to protect the PMDC motor and the 


power electronic converter. In this paper, limits were placed on the velocity reference by an amount that does 
not exceed the rated motor velocity, as well as limit, was placed on the reference voltage coming out of the PI 
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current by 5v, which represents the control voltage (vc) and is compared with a triangular voltage as 
explained [24]. 
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Figure 3. General structure 


3. TUNING METHOD 
3.1. Particle swarm optimization (PSO) 

Particle swarm algorithm (PSO) is proposed by Kennedy and Eberhart and it was modified to 
improve its performance by adding a new parameter called inertia weight [25]. It which bases on a swarm 
base for using an observation of social behavior of moving organisms Such as a gathering of fish or a flock of 
birds. 

The algorithm is related to the computational method that optimizes the problem Where flocks of 
birds aim to find eating behaviors and it uses repeated steps to reach the best solutions. This means that the 
candidate solutions are particles that move in the search space based on a specific formula above the particle 
position. Each particle's movement is affected by its local value, and its objective is to reach the best-known 
positions in the search-space by updating better positions found by other particles. 

The algorithm is begun by establishing the starting position and velocity vectors. And each iteration 
is tried to find the best value by evaluating position and velocity vectors. Every particle has variables and 
dimensions, and these variables are problems that need to be solved. If the problem consists of five different 
variables, the particles’ dimension should be chosen as five. This algorithm depends mainly on finding the 
position of each particle with the best local value, as well as finding the best general swarm position in each 
iteration. The position and speed are updated at each iteration based on (12) and (13) [21], {26]. 


Vege + 1) = W. Vij(t) T rıcı |Pbest; ;(t) = X,;(t)] F T2c2[Gbest; ;(t) = Xij (t)] (12) 
X,jt+1) = Vij(t+1) + X(t) (13) 


Where, i=particle index, t=iteration, Vaj (t+l)=velocity updated or new velocity, j=dimesion 
number, W=weighted inertia its value between [0-1], Vaj (t)=Current velocity, rı, r ~2=random coefficients 
there values between [0-1], Xaj (t)=current position, cı, co=acceleration coefficients there values between [0- 
2] and Xq, p (t+1)=position updated or new position. The particles in the PSO, position values are updated 
until the number of iterations is reached to specified value. The particles’ changing positions are illustrated in 
Figure 4 and Table 2 represents the parameters of the PSO that were used in this work. Figure 5 shows the 
flowchart of the particle swarm optimization algorithm. 


Table 2. PSO Parameters 


PSO parameters Value 
Iteration 10 
Swarm size 20 

No dimension 5 
Weighted inertia 0.9 
C1, C2 2, 1.5 
LB, UB 0, 300 
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point by the PSO [26] Figure 5. PSO flowchart [27] 


3.2. Butterfly optimization algorithm (BOA) 

In this paper, a new metaheuristic algorithm inspired by nature called the Butterfly Improvement 
Algorithm (BOA for short) has been proposed. This algorithm adopts the strategy of searching for food or 
mating between butterflies from a biological point of view. Butterflies use sense receptors to find the source 
of food by sensing the scent of the fragrance (nectar) [18], [28]. 

These receptors are spread over the body parts of the butterfly, such as the antennae and the legs, and 
they are nerve cells that are found on the surface of the butterfly’s body and are called chemical receptors. 
The butterfly generates a scent of different intensity according to its fitness, meaning that when the butterfly 
moves or moves from one place to another, the fragrance spreads across these distances, and other butterflies 
can sense it, and this is the method of communication between butterflies [18], [28]. 

The entire concept of sensing and odor processing is based on three important terms, a) the sensory 
modality, denoted by the symbol c, and its value between [0-1], b) stimulus intensity, symbolized by the 
symbol I, and c) the size of the stimulus or force on which the butterfly depends. It is known as power 
exponent and It is denoted by the symbol a and its value is between [0-1]. 

The natural phenomenon of butterflies is based on two important issues, which are the formulation of 
the fragrance function and the variation in the intensity of the fragrance. The fragrance function can be 
formulated as a function of physical density and according to (14). 


f = cl? (14) 


where f is the perceived magnitude of the fragrance, 1.e., how stronger the fragrance is perceived by other 


butterflies. c, I and a are explained by (14). 

When the butterfly can sense the smell emanating from another butterfly, it will move towards it, and 
this stage is known as the global search and can be represented by (15), and when the butterfly is unable to 
sense the smell emanating from any other butterfly, ıt will move randomly and this stage is called the local 
search and It can be represented by (16). 


Xi? = Xt + (r? * g*— Xi") «fj (15) 


Xt = Xi + (1? * Xt- X) *fi (16) 
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Where, X;‘*?=new solution vector, X;'=solution vector, r=rondom number (n) between [0,1], g=current best 
solution, t=iteration number, i=butterfly, X and Xx are jth and kth butterflies from the solution space, and 
fi=the perceived magnitude of the fragrance. 

The butterflies’ search for food or mating can occur on both local and global searches. The butterfly 
can switch between local search and global search, meaning that it has the decision to move towards the best 
butterfly or move randomly. The mode to switch between local and global search is called a switch 
Probability and denoted by the symbol p. Table 3 represents the parameters of the BOA that were used in this 
work. 

The butterfly optimization algorithm can be represented in the approximate flowchart shown in 
Figure 6. The design of any complex system requires the selection of certain criteria that give the best 
performance. These functions are known as performance indices. In this paper, ITAE was used and it can be 
represented by (17). 


ITAE = f; tle(t)|dt (17) 


Table 3. Parameters of the BOA 


Parameters Value 
Max iteration 5 
No search agents 20 | Set the valuesfor n , Lb , Ub, dim , ItrMax 
No demension 5 
Switch probability p 0.8 
POWPT ESD oneni a a Initialize the population of n butterflies X,‘ 
Sensory modality 0.01 ee | 
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LB, UB 0.300 l 
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Figure 6. BOA flowchart 
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4. RESULTS AND COMPARISON 
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Simulation was performed on the system using both the BOA and PSO methods. The system was 
tested in cases of loading and non-loading, as well as in the case of a single input or multiple inputs to ensure 
that the PMDC motor rotates a full rotation, 360 degrees or a half rotation. Table 4 represents the cascade PID 
parameters extracted by both the BOA and PSO methods, while Table 5 represents the performance 


parameter values obtained from these two methods. 


For the purpose of comparison between BOA and PSO: 

a. An overshoot was observed when using the PSO method and its value was 2.557%. This indicates a 
deviation of the position of 7.82 degrees from the reference position while no deviation of the position 
occurred when using the BOA method see Figure 7(a). 

b. An overshoot and velocity deviation of about 18rad/sec were observed when using the PSO, while no 
velocity deviation occurred when using the BOA method see Figure 8(a). 

c. When applying full load (17.6N.m) on the system at the fifth second, it was noticed that the system was 
not affected when using both methods, or was affected imperceptibly, and this indicates the durability of 
the cascade controller in rejecting the disturbance see Figure 7(b) and Figure 8(b) 

d. The BOA method reached the required values at the second iteration, while the PSO method needed 10 
iterations to reach the required values see Figures 9(a) and 9(b). 


Table 4. PID parameters values 

















PID parameters PSO BOA 
KD position 132.6568 77.3994 
KP velocity 102.3941 55.5583 
KI velocity 14.2140 0.524 
KP current 23.68556 74.3044 
KI current 0.1725 0.0236 
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Table 5. The values of the performance criteria for 
each tuning method 


Performance criteria PSO BOA 
Rise time (ms) 67.316 68.182 
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Figure 7 (a) Position control system at no load case (b) Close up load case position control at fifth second 
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Figure 8. (a) velocity control system at no load case (b) close up load case velocity control at fifth second 
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Figure 9. (a) BOA cost function vs iteration (b) PSO cost function vs iteration 


The simulation results showed that there was no position deviation or overshoot in velocity or 


position when using the parameters extracted by the BOA in addition to accurate tracking of the path. On the 
other hand, when using the PSO, a position deviation of 7.82 degrees and an overshoot of 2.557% were 
observed, in addition to an overshoot in the velocity. The simulation results and the comparison mentioned 
above showed that the BOA algorithm gave better results than the PSO algorithm. 


5. CONCLUSION 


In this paper, two optimization algorithms are used to extract the PID gains for the cascade controller 


of a PMDC motor. These algorithms are PSO and BOA. The simulation results in the case of using PSO 
showed a clear deviation of the position from the reference position about 7.82 degrees, in addition, an 
overshoot of 2.557%. On the other hand, no position deviation or overshoot occurred when using the BOA 
algorithm. Also, there is no overshoot in the velocity when using the BOA algorithm, which led to an 
accurate tracking of the path. As a result, the BOA algorithm is better than the PSO algorithm. 
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